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Vector mesons at finite temperature and QCD sum rules 
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Finite energy sum rules for vector and axial- vector currents are derived in a thermal medium to 
provide constraints for the spectral behavior of p and a\ mesons at nonvanishing temperature and 
hence to study the tendency toward chiral symmetry restoration. The parity-mixing ansatz for the p 
and a\ spectra, including finite widths, is investigated as a function of temperature. Characteristic 
differences between vector and axial- vector channels are discussed with regard to the implementation 
of the chiral-symmetry- breaking scale, 4n times the pion decay constant, in the sum rule approach. 

PACS numbers: 11.30.Qc, ll.30.Rd, 11.55.Hx, 14.40.Be 

Keywords: Vector meson, chiral symmetry, finite temperature and QCD sum rules 



1. INTRODUCTION 

The chiral SU(N f ) L x SU(N f ) R symmetry of the QCD 
Lagrangian with Nf (massless) quark flavors is known to 
be broken down spontaneously to SU(Nt)v in the non- 
perturbative QCD vacuum. With inclusion of explicit 
chiral symmetry breaking by the nonvanishing quark 
masses, this symmetry-breaking pattern is visibly man- 
ifest in the spectrum of the low-mass pseudoscalar me- 
son octet for Nf = 3. This octet figures as the mass- 
less Goldstone bosons of spontaneously broken chiral 
SU(3)l x SU(3)r symmetry in the limit of massless u, d, 
and s quarks. An equivalent but not directly observable 
manifestation of spontaneous chiral symmetry breaking 
(SxSB) is the appearance of a nonzero quark conden- 
sate. Its observable counterpart is the pion decay con- 
stant, /jr, which acts as an order parameter for S^SB and 
defines a chiral-symmetry-breaking scale, Awf^ ~ 1 GeV, 
briefly referred to as the chiral scale. Further evidence 
for S^SB is the absence of parity doublets in the lower 
mass area of the hadron spectrum, an important example 
being the large mass gap between the p (J p = 1~) and 
01 ( J p = 1 + ) mesons. Thus it is well established that at 
zero temperature, QCD is in the Nambu- Goldstone real- 
ization of chiral symmetry. Were it in the Wigner-Weyl 
realization with a trivial vacuum and parity doublets, the 
p and 01 would be degenerate. One expects this latter sit- 
uation to be reached as the temperature is raised toward 
a critical value, T c ~ 0.2 GeV, for the chiral transition. 
Around this point and beyond, the p and a\ are expected 
to melt into the quark-antiquark continuum. 

Lattice QCD results [H,[2] do indicate that the sponta- 
neously broken chiral symmetry is restored at high tem- 
perature. The spectral functions of chiral partners such 
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as the (isovector) vector and and axial- vector mesons are 
then supposed to become degenerate. In this context, 
concerning the in-medium behavior of vector mesons, 
measurements of dileptons produced in relativistic heavy- 
ion collisions have attracted great interest over the past 
decades (see, e.g., Ref. Q). Dileptons as well as photons 
are excellent probes carrying information about vector 
meson spectral functions at the high temperatures and 
densities reached in the collision zone. In particular, the 
p meson, as the lowest dipole excitation of the QCD vac- 
uum, has been in the focus of such investigations. 

The issue of in-medium hadronic properties persists as 
a basic theme ever since the suggestion of Brown and Rho 
(BR) Q that hadron masses should decrease systemati- 
cally when a hadronic medium undergoes the transition 
toward chiral symmetry restoration. The BR scaling hy- 
pothesis states that the dropping of hadron masses, with 
the exception of the pion mass, should go in parallel with 
the dropping of the pion decay constant f„ (here the one 
related to the time component of the axial current) in 
hot and dense hadronic matter. In general, a sufficient 
condition for chiral symmetry restoration is that the in- 
medium mass difference between the p and ai mesons 
tends to zero, although it is not mandatory that the 
masses vanish individually. In practice, however, these 
mass-shift scenarios are overshadowed by strong spectral- 
broadening effects. The p and a\ resonances have large 
widths already in the vacuum. Their interactions with 
hadrons in the medium increase their decay widths sig- 
nificantly through collision broadening. Distinguishing 
"mass shift" from "broadening" scenarios is obviously not 
a meaningful issue once the spectral functions become 
very broad and do not show clear resonance structures. 

In a situation like this, it nevertheless still makes sense 
to perform a sum-rule analysis of such broad spectral dis- 
tributions, focusing on their lowest moments and defining 
a mean mass through the first moment of the in-medium 
spectral function. Such an analysis is made possible by 
identifying the characteristic chiral-symmetry-breaking 
scale, Acsb ~ Airf^, with the continuum threshold sep- 
arating the low-energy-resonance region from the high- 
energy continuum. Consider the example of the vector 
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spectral function, to be specified in detail later: 

R v (s) = R r v es ( s ) 0(s v - s) + R c v ont Q(a - ay) , (1) 

where Ry s (s) denotes the resonance part. This low- 
energy part is separated by the scale sy from the high- 
energy continuum part Ry nt determined by perturba- 
tive QCD. In vacuum, the identification y / iy = 47r/ 7r 
is supported by current algebra relations and spectral 
sum rules, as we point out later. The extension to the 
case of finite baryon density at zero temperature has 
been studied in our previous work Q where it is demon- 
st rated t hat the aforementioned scale analysis works with 
\fsv\fi) = 47r/*(p) now interpreted in terms of the in- 
medium change of the chiral order parameter. 

Before turning to finite temperatures, it is useful to 
recall Weinberg's sum rules Q. These sum rules are en- 
tirely based on current algebra and establish rigorous re- 
lations for the difference between vector and axial- vector 
spectral functions: 



/>oo 

/ ds[R v (s)~ R A (s)] = 
Jo 

poo 

/ dss[R v (s) - R A {s)] = 
Jo 



(2) 
(3) 



Assume first the following schematic forms for the vector 
and axial- vector spectral functions: 



R v (s) = R r v es {s) Q(s v - s) + c 0(s - s v ) , 
R A (s) = R r A es (s) Q(s A - s) + c G(s - s A ) , 



(4) 



where Co multiplying the continuum parts of the spectral 
functions is determined by perturbative QCD (co = 3/2 
in leading order for the p and eti channels). The equal- 
ity of these vector and axial-vector continuum pieces is 
indicative of chiral symmetry restoration at high energy 
beyond the scales sy,A- For the resonant parts, choose for 
simplicity a zero- width ansatz (as realized in the large- N c 
limit) : 



R r y es ( S ) = 12tt 2 f v m v 8(s - m v ) , 

RT(s) = 127T 2 [f 2 8(s - ml) + j\m\ S(s - m|)] 



(5) 



where fy/ A are dimensionless vector and axial- vector 
couplings respectively. The a\ spectral function includes 
a pion pole term with its residue determined by the pion 
decay constant /„.. In the chiral limit (m^ — > 0), and 
assuming equal continuum thresholds for the vector and 
axial- vector channels, sy = sa, the Weinberg sum rules 
(12I3D read as follows: 



t2 2 j-2 2 t2 
fv m V ~ JA m A = L , 
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fv m v 



flm\ = 



(6) 



These equations, together with the KSRF relation 0,Q, 
f v m v = 2/ 2 , (7) 



imply that the vacuum a\ and p masses satisfy the well- 
known relation 



2m v . 



(8) 



The Weinberg sum rules are based on the observation 
that chiral symmetry is restored in its Wigner-Weyl real- 
ization asymptotically, at high-energy scales, where QCD 
is perturbative. The duality between the resonant and 
asymptotic parts of the spectral functions derived from 
current correlators is one of the basic ideas of QCD sum 
rules, to which we now proceed. The present study con- 
structs finite energy sum rules (FESR) for p and <xi me- 
son at nonzero temperature and vanishing baryon chem- 
ical potential. We show how finite temperature FESRs 
can provide constraints for the pattern and trend toward 
chiral symmetry restoration. Realistic vacuum spec- 
tral functions are employed as input in comparison with 
schematic 8 function spectra, and the p-a\ parity-mixing 
scenario is used to describe temperature-dependent ef- 
fects. 



2. FINITE ENERGY SUM RULES 

Following these preparations, the starting point is now 
the time-ordered current correlation function 



W v {q)=i / d A xe^ x (Tf{x)f{Q)h 



(9) 



with the vector current j v {x) = \ (yrf^u — d/y^d) and 
the axial- vector current = i (uj^-f^u — dj^j^d) car- 
rying the quantum numbers of p and a\ meson, respec- 
tively. The bracket (C)t indicates the thermal expecta- 
tion value of an operator 0, 



(Oh 



trO exp (-H/T) 
trexp (-H/T) ' 



(10) 



where H is the Hamiltonian. 

In vacuum, the tensor correlation function (j9]) can be 
related to a single invariant correlator, H(q 2 ) = ^g^H^ . 
In a medium, the longitudinal and transverse parts of the 
correlator differ as a consequence of broken Lorentz in- 
variance. In the rest frame of the medium and for the case 
in which the mesons have vanishing three-momentum 
q = 0, however, longitudinal and transverse correlation 
functions coincide and again are given as a single func- 
tion H(uj,q = 0). This correlator is written in the form 
of a twice-subtracted dispersion relation: 



n( 9 2 ) = n(o) + n'(o)q 2 + ^ J di 



Imh(s) 



s 2 (s — q 2 — ie) 



(11) 



On the other hand, the operator product expansion 
(OPE) is used to represent the correlator at large space- 
like momentum, q 2 = — Q 2 < 0: 



12^ 2 nV = -Q 2 ) = -c Q 2 ln ( ^j+ Cl+ ^ + ^- 



(12) 
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with the coefficients 



co = 2 i 1 + M , 
9 

7T 2 Ct 

c 2 = — (— G 2 ) ± 67T 2 (m u (uu) T + m d (dd) T ) 



(13) 



Here, 5^ in Co denotes the radiative corrections in per- 
turbative QCD. The explicit form of 5 N up to order a^(s) 
can be found in Ref. j5|- The difference between vector 
and axial-vector channels results to this order from the 
sign of the quark condensate term in c-2- In the chiral 
limit (m Mj d — > 0), c\ and the second term of Ci van- 
ish. The first appearance of a difference between p and 
ai spectral functions is then in the term proportional to 
C3 in the expansion (|12[) . This term involves four-quark 
condensates, 



q—u,d,s 



(14) 



that are subject to large uncertainties. To evaluate these 
condensates, a factorization approximation is frequently 
used, assuming that intermediate states are saturated by 
the QCD ground state: 
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m^ 5 X a q) 2 } = -((q^X-qY) = - K (qq) 



(15) 



with k introduced to parametrize the deviation from ex- 
act factorization (k = 1). 

However, the factorization approximation for the four- 
quark condensates is not sufficiently accurate for any 
quantitative considerations This basic uncertainty 
can be avoided by arranging sum rules in terms of mo- 
ments of the spectral function and restricting the analysis 
to the lowest two moments: 



/ dsR{s,T) 
Jo 

so CO c 



»o(T) 



ds s R(s, T) 



ci - 12tt 2 n(0) , (16) 
^ co -c 2 (T). (17) 



The dimensionless, temperature-dependent spectral 
function R(s, T) stands generically for the vector or axial 
channel and is defined as 

R(s,T) = -— Imn(s,T) . 
s 

Its high-energy continuum is separated from the low- 
energy part by the scale so(T) that is expected to shift 
downward with increasing temperature. Again, this sq 
stands for the continuum threshold scale in either vector 
or axial-vector channel. These scales, denoted in the fol- 
lowing by Sv(T) and sa(T), are generally different and 



to be determined by the detailed sum-rule analysis. The 
last term on the right-hand side of Eq. ([T6"]) vanishes in 
the vector channel and represents the pion pole contri- 
bution in the axial-vector channel, with n(0) = / 2 at 
T = 0. 

In the asymptotic region (Q 2 — > 00) where the OPE 
is valid, all nonperturbative scales appear in the form of 
power corrections to the perturbative calculations. These 
nonperturbative contributions are separated into the re- 
spective condensates. At low temperatures, it can be 
assumed that, apart from the scales sv,a{T), the T de- 
pendence is only in the condensates. 



3. EVALUATION OF (0) T 

The evaluation of the T-dependent condensates follows 
the method employed previously in Ref. [j| . In Eq. (|TU1) , 
the vacuum state and the lowest excitations of the hadron 
gas, namely pions, are taken into account as eigenstates 
of the Hamiltonian to calculate the thermal expectation 
values at low temperatures: 



(0) T = (0) Q +it J 



d 3 p 



2E{2ir) 



(<K a (p)\0\ir a (p))n 



B j 



where n 



D 



,E/T 



1 



(18) 

denotes Bose-Einstein distri- 



butions of thermal pions with E 



+ p 2 . The pion 



matrix element in Eq. 
pion limit : 

(7T a (p)\OK(p))=- 



181) can be evaluated in the soft 



f 2 

J 7T 



[eg, [eg, o}} 10) 



(19) 



where Q5 is the axial charge operator. Equations (fl~8|) 
and (fT9|). when applied to the scalar quark operator, qq, 
give the leading order T dependence of the chiral con- 
densate, 



with 



(qq)T = (qq)o 



B^x) 



1 



T 2 



T 



dy 



ev - 1 



(20) 



(21) 



Equation ([2H)l is the well-known result derived from chiral 
effective field theory 

The application of Eq. ([T5]) to the gluon condensate 
makes use of the QCD trace anomaly, 



6 M = — 
8 



11 - -N f 

3 } 



-G 2 + ^m 9 g g , (22) 

9 



to calculate the pionic matrix element of the relevant 
gluon operator [13. For Nf = 3, the resulting tempera- 
ture dependent gluon condensate becomes 



—G 2 ) 

7T IT 



-G" 



1 



-mlT 2 B 1 



run 
T 



(23) 
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in which the second term of the right-hand side gives 
a numerically minor contribution to the sum rule. In 
the actual calculation, the Gell-Mann— Oakes— Renner 
(GOR) relation is used for evaluating the vacuum quark 
condensate, m„ (qq) = — (O.llGeV) 4 , while the charmo- 
nium sum rules [l4( constrain the vacuum gluon conden- 
sate as (^G 2 ) = 0.005 ± 0.004 GcV 4 , with large uncer- 
tainty. In practice, this uncertainty is not prohibitive 
since the overall correction induced by C2 in Eq. (fT7|) 
is small compared to the dominant term proportional to 
4(T). 

Apart from the in-medium modifications of the scalar 
condensates, new spin-dependent operators appear in the 
OPE because of broken Lorentz invariance in the heat 
bath. Such operators are classified by their canonical di- 
mension and twist (r = dimension — spin) . The C2 (T) in 
the sum rule for the first moment now includes the con- 
tribution from the twist-2 operator ST q^/^D^q, where 
the symbol ST makes the operator symmetric and trace- 
less with respect to its Lorentz indices. The pion matrix 
element of this operator in Eq. (|T9"|) is evaluated using its 
relation to the quark distribution function in the pion Q . 
In the actual calculation, this contribution to C2(T) is 
about three times smaller than the leading T-dependent 
correction to the quark condensate in Eq. (|20[) . 



4. PHENOMENOLOGY AT LOW T 

The mixing of vector and axial- vector correlation func- 
tions at finite temperature is introduced following Ref. 
(l5| and translates correspondingly to the spectral func- 
tions: 

R v (s,T) = [l-e(T)] R v (s,0) + e(T)R A (s,0) , 
R A (s,T) = e(T)R v (s,0) + [l-e(T)] R A (s, 0). (24) 

The mixing parameter e is given by the thermal pion loop 
integral 



e(T) 



2 

72 

■/ 7T 



V 



1 



(2tt) 3 E(e E / T - 1) 



(25) 



In the chiral limit (m n — > 0) it reduces to e = T 2 /(6f 2 ). 
At the temperature Td where e ~ 0.5, the vector and 
axial- vector spectral functions are maximally mixed and 
become degenerate, Ry{s,T) = R A {s,T). The tempera- 
ture at which this degeneracy takes place, Td = %/3/tt — 
151 MeV (with f n = 87 MeV in the chiral limit), is not 
far from the characteristic temperature commonly asso- 
ciated with the chiral crossover transition (T c = 170 — 190 
MeV) P, Q| . Of course, as one approaches this tempera- 
ture range, terms of higher order in T 2 become important 
and need to be taken into account. 

It was pointed out in Ref. [HI that the vector meson 
mass changes weakly at low temperature (i.e., at order 
T 2 ). We demonstrate this feature next by performing the 
FESRs introduced in the previous sections. 



A. Spectral functions with zero width 

It is instructive to test the FESR for the vector and 
axial-vector spectral functions with zero width, the lim- 
iting situation realized in the large N c limit of QCD: 



R v (s,T = 0) 



12tt 2 f 2 m 2 v 5(s 



R A (s,T = 0) = 12Tr 2 f 2 



5{s 



4) 



(26) 



The vector and axial- vector couplings, fv/ A , are deter- 
mined by 



f 2 m 2 v = 2fl f A m\~fl 



s A - s v 



(27) 



The first part is just the KSRF relation; the second one 
is derived by applying the KSRF relation to the first 
Weinberg sum rule in Eq. © and adding the correction 
coming from the difference between s A and sy. Then the 
lowest two moments of the vector sum rule are: 



12tt 2 f v m v [l-e(T)} = 



12k 2 f 2 m v [I- e{T)] = 



-s v (T)-l2n 2 f 2 e(T), 



^s 2 v (T)-c 2 (T) , 



(28) 
(29) 



where the a\ meson contribution has been omitted be- 
cause in the actual calculation, y/sy is found to be lo- 



cated at sy < m 2 A so that the a\ pole does not contribute 
to the integrals of the spectral moments. 

With this setup, the continuum threshold in the vector 
channel is evaluated to be ^/~sy — LI GeV at T = from 
the sum rules, using f v ~ 87 MeV in the chiral limit and 
ray = 770 MeV. Note that yfsy turns out to be perfectly 

close to the chiral scale, An f^. 

The finite temperature behavior of yj sv(T) extracted 
from the sum rules to leading order in e is shown in Fig.[T] 
The sum rules for the lowest two moments give consistent 
results at low temperature. However, although the iden- 
tification y / sy = 47r/ T emerges naturally from Eq. (|28p at 
T = 0, the temperature dependence of y/sy evolves more 
slowly than the canonical f n (T) = [1 - e(T)/2] / W (T = 0) 
found in chiral perturbation theory (Tp| . 

Defining an average mass fh by the normalized first 
moment of the spectral function, 



m 



J ds sR 
JdsR 



(30) 



we confirm that the p meson mass remains unchanged 
for T^0 and stays at the vacuum pole position, rhy = 
770 MeV, to leading order in e. This is consistent with 
the statement f]~5l. Il~6{| that the vector meson pole position 
remains unchanged at order T 2 . 

In contrast to the p meson case in which the a\ pole 
does not contribute to the integrals of the spectral mo- 
ments (since sy < m 2 ^), the left-hand side of the sum 
rules in the axial-vector channel receives contributions 
from both p and a\ poles. By using Eq. (|27|) the sum 
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FIG. 1: (Color online) Continuum threshold yfsv in the vec- 
tor meson channel as a function of T obtained from the sum 
rules for the 0th (black solid line) and 1st (red dashed line) 
moment of the spectral function. The T dependence of chiral 
scale 47r/^(T) (blue dotted line) is also displayed. 
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FIG. 2: (Color online) Axial-vector continuum threshold ^/sA 
as a function of T obtained from the sum rules for the 0th 
(black solid line) and 1st (red dashed line) moment of the a\ 
spectral function. In contrast to the p meson, the ai-meson 
mass (black dotted line) decreases with rising temperature. 



rule for the first moment of the a\ spectral function be- 
comes 

127T 2 [f\ m\(l - e) + f v m 4 y e] = \s 2 A - c 2 . (31) 

With m A = 1.26 GeV, this gives ^/sa ^1.54 GcV at 
T = and thus f\ ~ 0.014, which agrees with the em- 
pirical values of f\ in Ref. [l7j. As seen in Fig. [2J the 
average axial-vector mass fh,A defined by Eq. ([30]) de- 
creases in parallel with the continuum threshold y/sA as 
the temperature rises. One observes that the cti mass 
fa a — 1-09 GeV does not yet approach the p mass 
my = 770 MeV at temperatures around Td — 151 MeV, 
where e = 0.5, suggesting that higher order effects in e 
become important as one gets closer to the critical region 
for chiral restoration. 



B. Realistic spectral functions with finite width 

In this section we study realistic empirical spectral 
functions with finite widths as displayed in Fig. [3] The 
continuum is still described by a step function. In pre- 
vious work Q, we have shown that, both in the vac- 
uum and at finite baryon density, the detailed choice for 
modeling the continuum threshold is not decisive for the 
outcome of the sum-rule analysis. In practice, the con- 
tinuum threshold can be given an improved description 
using a ramping function, resulting in the blue solid curve 




0.0 0.5 1.0 1.5 2.0 2.5 3.0 



s [GeV 2 ] 

FIG. 3: (Color online) Vector (black curve) and axial-vector 
(red curve) spectral distributions in vacuum, compared to 
e + e~ — > nn data with n even fl8l. fTglj and data from hadronic 
r decays (20I. [2l|. Here sv and sa stand for the continuum 
thresholds in vector and axial-vector channels, respectively. 
The ramping function (blue dashed line) shows an example of 
smooth threshold modeling (see text and appendix). 

in Fig. [3l which reproduces the experimental data very 
well. We show that for the finite-temperature sum rules 
considered here, the results are also consistently stable, 
independent of the detailed threshold modeling, as long 
as the slope of the ramping into the continuum is chosen 
sufficiently large (see Appendix [A"|) . 
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FIG. 4: (Color online) The p and a\ meson spectra at finite 
temperature. 
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Figure^] exhibits the temperature-dependent p and a\ 
spectra generated by Eq. ([24)) with physical values of the 
pion mass and decay constant, = 139.6 MeV and 
fix = 92.4 MeV. The pole position of the vector me- 
son resonance stays at its vacuum mass. However, when 
examining the average mass defined by Eq. (|30[) , 



f4(T)-c 2 (T) 



c s v (T) + Cl - 12^ 2 /2 e(T) 



(32) 



the broad widths of the p and a± spectra affect the p 
mass distribution. From Table U] it is apparent that the 



T (MeV) 


rh v (MeV) 


V^V (GeV) 





791 ± 2 


1.139 ±0.007 


50 


792 ±2 


1.138 ±0.007 


100 


800 ±4 


1.133 ±0.008 


120 


807 ±5 


1.128 ±0.009 



TABLE I: Average vector mass and continuum threshold at 
various T from Eqs. (|32[1 and (|17[) . The errors include un- 
certainties in the value of a s (entering the NLO perturbative 
QCD corrections) and of the gluon condensate. 



average vector mass tends to slightly move upward with 
temperature, picking up weight from the mixing with the 
Oi, while y/~sv moves downward. When using the tightly 
constrained vacuum spectral functions as input at T = 0, 
it turns out that there can be a small (<5%) mismatch 
between the left-hand and the right-hand sides of the 
sum rules. The error assignments in the values of y / sy" 
from the sum rule for the first spectral moment take this 
small uncertainty into account in Table[l] In practice, the 
errors in Table |T] come primarily from the uncertainties 
of the gluon vacuum condensate and the running strong 
coupling a s (s). 

As the temperature increases the continuum onset 
scale sy (T) shows evidentally a tendency to decrease. At 
temperatures T > 140 MeV, however, solutions for ^/sy 
cease to exist, and hence the average mass at that high 
T cannot be determined via Eq. (|32[) . This is partly a 
consequence of the restrictions imposed by the treatment 
to leading order in e(T). Even more so, it is related to 
the fact that, with realistic spectral functions, the broad 
a\ distribution does not permit a separation of scales be- 
tween low- and high-energy parts of the spectrum. The 
increased mixing of the ai into the p spectrum at high 
temperature enhances this effect. Thus, with realistic 
spectral distributions, the FESR approach has its ap- 
plicability limited to the range well below the "melting" 
temperature of the resonances. 

As mentioned in the introduction, the difference be- 
tween vector and axial- vector correlators serves as an or- 
der parameter of spontaneous chiral symmetry breaking 
and restoration. Direct substraction of the FESRs in the 



axial-vector and vector channels gives 

s V (T) rSA (T) 

dsR v (s,T) - / dsR A (s,T) 
o Jo 

= c [s v (T) - s A (T)] ± 12tt 2 /' [1 - 2e(T)] , 

sv{T) r s A (T) 

ds s R v (s, T) — dss R A (s, T) 



CO ln2 

2 



[4(T)-4(T)] 



(33) 



Keeping the leading order mixing ansatz (|24[) for the 
spectral functions, one finds ~ 1.12 GeV and ^/sa — 
1.49 GeV at T = from Eqs. ([53]). Note that „/sy~ is 
again close to the chiral scale 47r/ ff . It is important to 
realize that sa, as it results from the consistency con- 
ditions of the sum-rule analysis, is different from (i.e., 
significantly larger than) sy. In fact, the a\ mass itself is 
comparable to 47r/ T so that the continuum threshold in 
the axial- vector channel must necessarily be located at a 
higher scale than sy once the large a\ width is taken into 
account. Assuming sa — sy = sq as in the schematic 
Weinberg sum rules ©, one finds that the leading or- 
der temperature dependence, 1 — 2e, just drops out in 
Eq. (|3"3"j) and yields a trivial result, 



/ °ds [i?v(s,0)-i? A ( S ,0)] = 127r 2 /' , 
Jo 

/ ds s [R v {s, 0) - R A (s, 0)}=Q , 
Jo 



(34) 



independent of temperature. 

5. HIGHER ORDER CORRECTIONS 

A. T 4 corrections 

Effects at order T 4 have been treated systematically in 
Ref . [l|| . It is shown there that these corrections amount 
to replacing the mixing parameter e by e — > e(l — e/2). 
The temperature T ( i at which the p and a\ spectral dis- 
tributions become degenerate is now shifted upward to 
Td = 

The T correction incorporates interacting pions with 
nonvanishing momentum in the heat bath, which con- 
tribute to the Lorentz noninvariant part of the current 
correlation function. According to the Ref. [l||, this con- 
tribution can be expressed in the frame with q = as 
follows: 

Tly(q 2 = -Q 2 ,T) = Tly(-Q 2 ,T = 0) 

- e (l - |) [Uy(-Q 2 , T = 0) - n A (-Q 2 , T = 0)] 



4^ 2 M 2 4 
15 Q 2 



(35) 
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where M2 is the first moment of quark distributions in 
the pion: 



6. SPECTRAL FUNCTIONS FROM AN 
EFFECTIVE FIELD THEORY 



Mo 



dxx [v(x) + 2s(x)] 



(36) 



with valence and sea quark distributions, v(x) and s(x), 
resp ectively. We use the value M2 ~ 0.12 as discussed in 

By transferring Eq. (|3"5|) to the FESR, it turns out 
that Eq. (fT7|) is modified by the new term of order T A as 
follows: 



:v(T) , , x 16tt 4 M 2 , 
dssR v {s,T) ^~=-T 4 



|4(T)-c 2 (T), 



(37) 

while the denominator of Eq. (|32p receives a small cor- 
rection from the replacement e — > e(l — e/2). 

On the OPE side, the T dependence of the quark con- 
densates is also improved up to order T 4 . Chiral pertur- 
bation theory gives [ll| 



1 ~ \ e ~ J> £2 



(38) 



However, this correction is numerically small in the ac- 
tual calculation. 

The substitution of Eq. (|3"7| in the numerator of 
Eq. (I32p introduces a relatively small negative mass shift, 



5mv 



5/ n " (T) dsR v (s,T) 



(39) 



consistent with the findings of Ref. [l6| , which compen- 
sates for the increase of the vector average mass by the 
finite width effect. For instance the average mass at 
T = 100 MeV is now obtained as m v ~ 798 ± 4 MeV, 
slightly less than in Table U but only marginally different. 



So far we have used empirical input for the spectral dis- 
tributions at zero temperature. In this section, we illus- 
trate how the sum rule works for the example of a chiral 
effective field theory in which the vector meson spectral 
distribution can be explicitly calculated and tested with 
respect to its sum-rule consistency. The example chosen 
is the chiral Lagrangian based on generalized hidden lo- 
cal symmetry (GHLS) [2^|, which explicitly includes the 
axial-vector meson in addition to pion and vector me- 
son. Details of the formalism at one loop are found in 
Ref. [13. 

The parity mixing in the p meson spectrum at finite 
temprature is generated from a process in which 7r and a\ 
mesons circulate in a loop attached to the p meson, with 
the pion coming from the heat bath. In this approach 
one does not have to rely on equations such as Eqs. (pM)) 
or J33J. 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 1.4 
s [GeV 2 ] 

FIG. 5: The p meson resonance at finite temperature calcu- 
lated in the GHLS model. 



B. Massive states 

For higher temperatures the contributions from mas- 
sive excitations such as K and r] in Eq. (JTUJ) need to 
be considered. In the T dependence of the quark and 
gluon condensates, these contributions are included as in 
Ref. 0: 




In practice, the T dependence of the gluon condensate 
is just a few percent of its vacuum value and negligible. 
These numerically minor corrections in the OPE are of 
only little significance to the sum-rule analysis. 



Figure [5] shows the vector spectral function of the 
GHLS approach at several temperatures. The spectrum 
exhibits the usual resonant peak, the width of which is 
governed by the p — > 2-k decay. The p-a\ mixing in- 
volves mechanisms in which the energy of the timelike 
p meson splits into two branches corresponding to the 
processes p + 7r — > a\ and p — > a\ + tt, with thresholds 
■s/s = m ai — run and yfs = m ai + m v . This produces the 
threshold effects seen as a shoulder at = m ai — m w 
(and a bump above ^/s = m ai + m w , not seen in the 
figure). The height of the p spectrum gets reduced with 
increasing temperature, whereas the ai-meson contribu- 
tion is enhanced via the mixing effect. 

The average vector mass defined in Eq. (f3TJ)) and the 
continuum threshold in the vector channel as a function 
of temperature are summarized in Table ITT] One ob- 
serves again that my stays unchanged from its vacuum 
value over the whole temperature range T < 140 MeV. 
The continuum threshold scale ^/sy shows a systematic 
decrease toward higher temperature, but at a rate con- 
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T (MeV) 


fhv (MeV) 


yfiV (GeV) 





787 ±2 


1.126 ±0.007 


40 


787 ±2 


1.126 ±0.007 


80 


787 ±2 


1.122 ±0.007 


120 


786 ±3 


1.111 ± 0.008 


140 


786 ±3 


1.102 ±0.008 



TABLE II: The average p-meson mass and the continuum 
threshold i/sv at various temperatures, following from the 
FESR consistency test of the spectral function calculated in 
the GHLS approach [H|. 

siderably smaller than the expected behavior of the chiral 
order parameter, f n (T). 

7. CONCLUSIONS 

In this work we have constructed FESRs in order to 
study the behavior of p and a± mesons as well as their 
mixing at finite temperature, with the aim of exploring 
the pattern of chiral symmetry restoration. The sum 
rules for the lowest two spectral moments of vector and 
axial-vector spectral functions involve only the leading 
QCD condensates as corrections. With inclusion of per- 
turbative QCD terms up to order a,, these sum rules 
permit a reliable quantitative analysis, unaffected by the 
large uncertainties from condensates of higher dimension 
such as the four-quark condensates. 

The leading temperature corrections involve thermal 
pion loops. To order T 2 , the temperature dependence is 
generated entirely by the p — a\ mixing effect caused by 
the nonvanishing pir a\ coupling in the presence of a pio- 
nic heat bath, as pointed out previously by Dey, Eletsky, 
and Ioffe (l^. To this order, there is no shift of the pole 
mass in the thermal vector meson propagator. Next-to- 
leading effects of order T 4 involve ttit intermediate states 
and can lead to moderate mass shifts. 

In the large- N c limit the p and a\ widths vanish. In this 
limit, schematic distributions with 8 function resonances 
and a step-function parametrization of the high-energy 
continuum, when inserted in the FESR analysis, repro- 
duce the well-known current algebra and chiral sum rules. 
As an interesting feature, one finds that at zero tempera- 
ture the continuum threshold sy in the vector channel is 
identified with the scale characteristic of spontaneous chi- 
ral symmetry breaking: y / iy = 4n f v . Above this scale, 
chiral symmetry is restored in its Wigner-Weyl realiza- 
tion. In the resonance region below this scale, the sym- 
metry is in the spontaneously broken Nambu-Goldstone 
realization. 

When realistic spectral functions with large widths are 
implemented, this "clean" separation of scales is no longer 
rigorously maintained, but the FESR analysis is still use- 
ful, with the pole mass in the vector correlator now re- 
placed by the normalized first moment of the spectral 
distribution. Within its range of applicability (up to 



temperatures of about 140 MeV), the sum-rule analy- 
sis consistently shows an almost constant behavior of 
this average vector meson spectral mass. The primary 
temperature dependence of the spectral function comes 
from p-a\ mixing in the thermal pionic medium. The a\ 
mass, again identified with the normalized first moment 
of the axial-vector spectral distribution, decreases with 
rising temperature. This indicates the expected tendency 
of the p and a\ spectra becoming identical (degenerate) 
when chiral symmetry is restored. However, the present 
analysis does not support the BR scaling hypothesis of a 
dropping p meson mass, at least not up to temperatures 
T ~ 140 MeV, which are not far from the chiral crossover 
transition temperature, 170 — 190 MeV. 

The continuum threshold scale sy in the vector chan- 
nel (even when smoothed by a ramping function) system- 
atically moves downward in energy as the temperature 
increases. This feature is observed in all cases studied. 
However, while the identification of y/sy with the chiral 
scale 47r f„ emerges naturally at T = 0, the downward 
evolution with temperature of y/sv is significantly slower 
than that of / w (T) deduced from chiral perturbation the- 
ory. 

ACKNOWLEDGEMENTS 

This work has been supported in part by BMBF, GSI, 
and the DFG Cluster of Excellence Origin and Structure 
of the Universe. Two of us (WW and CS) thank the or- 
ganizers of the program 'New Frontiers in QCD' for their 
kind hospitality at the Yukawa Institute of Theoretical 
Physics, Kyoto, where this article was finalized. YK is 
grateful to S. H. Lee and S. i. Nam for useful discus- 
sions. CS acknowledges parital support by the Hessian 
LOEWE initiative through the Helmholtz International 
Center for FAIR (HIC for FAIR). 

Appendix A: Continuum threshold modeling 

Here we test the reliability of the continuum threshold 
parametrization by a schematic step function. Such a 
test can be performed by replacing the step function with 
a ramp function to yield a smooth transition between 
resonance and continuum regions, as follows: 

R(s) = R res (s) 6(s 2 - s) + R c (s) W(s) , (Al) 

where the weight function, W(s), is defined as 

{0 for x < Si 
X ~ Sl for si < x < s 2 (A2) 
1 for x > S2 • 

The step function behavior is recovered for W(x) in the 
limit si — > S2- 
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1.155 




Slope of W(s) [GeV" 



FIG. 6: Dependence of yfsv [determined from Eqs. (|A3|I - 
(| A5[) ] on the slope (s2 — si) -1 of the ramp function W(s) 
describing the onset of the continuum for the example of T = 
100 MeV. The grey band indicates the uncertainty range of 
the result obtained with step-function parametrization of the 
continuum. 



By using the function W(s), the modified sum rules for 
the lowest two moments of the spectrum R(s) become 



f 

Jo 



dsR res (s) = s 2 ( c + ^£ ) + ci - 127r 2 n(0) 



[CO -RresM] [ 2 dsW(s) 
J si 



(A3) 



f S2 s 2 ( 3 

J ds sR res {s) = — ( c + -Ei 



C2 

S2 



[CO - Rres(s 2 )] f ' ds sW(s) . 
J si 



(A4) 



Sets of intervals [s±, 82] are then determined so as to sat- 
isfy both sum rules (|A31 IA4|) . and the scale sy, defined 
by 



sv 



si + s 2 



(A5) 



is now introduced to characterize the continuum thresh- 
old. In any temperature region where the sum rules are 
valid the present analysis does not depend sensitively on 
details of the threshold modeling. For an example of the 
p meson channel, it is demonstrated in Fig. [6] that the re- 
sulting sfsv at finite temperature is stable with respect 
to variations in the slope (s2~si) 1 of the ramp function 
W(s). In this test the uncertainties of a s (Q 2 ) and of the 
gluon condensate have been disregarded for simplicity. 



[1] M. Cheng et al, Phys. Rev. D 74, 054507 (2006); Phys. 
Rev. D 77, 014511 (2008); Phys. Rev. D 81, 054510 
(2010). 

[2] Y. Aoki, Z. Fodor, S. D. Katz and K. K. Szabo, Phys. 
Lett. B 643, 46 (2006); Y. Aoki et al., JHEP 0906, 088 
(2009). 

[3] R. Amaldi et al. (NA60 collaboration), Eur. Phys. J. C 
49 (2007) 235. 

[4] G. E. Brown and M. Rho, Phys. Rev. Lett. 66, 2720 
(1991). 

[5] Y. Kwon, M. Procura and W. Weise, Phys. Rev. C 78, 

055203 (2008). 
[6] S. Weinberg, Phys. Rev. Lett. 18, 507 (1967). 
[7] K. Kawarabayashi and M. Suzuki, Phys. Rev. Lett. 16 

(1966) 255. 

[8] Riazuddin and Fayyazuddin, Phys. Rev. 147 (1966) 
1071. 

[9] T. Hatsuda, Y. Koike and S. H. Lee, Nucl. Phys. B 394 
(1993) 221. 

[10] J. Gasser and H. Leutwyler, Phys. Lett. B 184, 83 (1987). 
[11] P. Gerber and H. Leutwyler, Nucl. Phys. B 321, 387 
(1989). 

[12] N. Kaiser, Phys. Rev. C 59, 2945 (1999). 

[13] T. Hatsuda and T. Kunihiro, Nucl. Phys. B 387, 715 



(1992). 

[14] B. L. Ioffe, Prog. Part. Nucl. Phys. 56, 232 (2006). 
[15] M. Dey, V. L. Eletsky and B. L. Ioffe, Phys. Lett. B 252, 
620 (1990). 

[16] V. L. Eletsky and B. L. Ioffe, Phys. Rev. D 51, 2371 
(1995). 

[17] J. J. Sakurai, K. Schilcher and M. D. Tran, Phys. Lett. 

B 102, 55 (1981). 
[18] A. Aloisio et al. [KLOE Collaboration], Phys. Lett. B 

606, 12 (2005). 
[19] S. I. Dolinsky et al, Phys. Rept. 202, 99 (1991). 
[20] R. Barate et al. [ALEPH Collaboration], Eur. Phys. J. C 

4, 409 (1998). 

[21] K. Ackerstaff et al. [OPAL Collaboration], Eur. Phys. J. 

C 7, 571 (1999). 
[22] M. Bando, T. Kugo and K. Yamawaki, Nucl. Phys. B 

259, 493 (1985); Phys. Rept. 164, 217 (1988); M. Bando, 

T. Fujiwara and K. Yamawaki, Prog. Theor. Phys. 79, 

1140 (1988); N. Kaiser and U. G. Meissner, Nucl. Phys. 

A 519, 671 (1990). 
[23] M. Harada and C. Sasaki, Phys. Rev. D 73, 036001 

(2006); M. Harada, C. Sasaki and W. Weise, Phys. Rev. 

D 78, 114003 (2008). 



9 



